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Metal-insulator transitions and diflerent ground-state phases in quasi-one-dimensional materials, 
' (_Riii2-DCNQI)2M (_Ri=_R2=CH3, I and M=Ag, Cu), are studied with a renormalization-group 

' method. We use one-dimensional continuum models with backward scatterings, umklapp processes 

and couplings with 2kF and ikp phonons (not static lattice distortion). We take a quarter-filled 
band for M=Ag and a sixth-filled band coupled with a third-filled band for M=Cu. Depending 
on electron-electron and electron-phonon coupling strengths, the ground-state phase becomes a 
Tomonaga-Luttinger liquid or a state with a gap(s). For M=Ag, there appear a spin-gap state with 
a dominant 2kF charge-density-wave correlation, a Mott insulator with a dominant Akp charge- 
density-wave correlation, or a spin-Peierls state with different magnitudes of spin and charge gaps. 
Three-dimensionality is taken into account by cutting off the logarithmic singularity in either the 
particle-particle channel or the particle-hole channel. The difference between the ground-state phase 
^ , of the f?i=7?2=CH3 salt (spin-Peierls state) and that of the Ri=R2=l salt (antiferromagnetic state) 

is qualitatively explained by a difference in the cutoff energy in the particle-particle channel. For 
M=Cu, there appear a Mott insulator with a charge density wave of period 3 and a Peierls insulator 
^ . with a charge density wave of period 6. The conditions for the experimentally observed, Mott 



oo 



insulator phase are strong correlation in the sixth-filled band, moderate electron-phonon couplings, 
and finite electron-4fcF phonon coupling. Resistance is calculated as a function of temperature with 
a memory-function approximation in both cases above. It qualitatively reproduces the differences 
' among the M=Ag and M=Cu cases as well as the -Ri=7?2=CH3 and Ri=R2=l cases. 

c ' 
o 

■ I. INTRODUCTION 

> 

m 

oo ' Quasi-one-dimensional materials, (i?ii?2-DCNQI)2M (M=Ag, Cu), show various phases. For M=:Cu, the hybridiza- 
T-H tion of DCNQI vr orbitals and Cu d orbitals causes intriguing transport 0] and magnetic |^ properties. The metal- 
insulator transition is accompanied by the formation of a charge density wave (CDW) of period 3. |^,^ The low- 
temperature phase is a Mott insulator, which is caused by strong correlation of d electrons in cooperation with 
electron-lattice coupling and cannot be described by a band picture. In fact, the spin susceptibility in the insulating 
state is enhanced over the temperature-independent Pauli-like paramagnetism seen in the metallic state. It has called 

■ forth theoretical studies, most of which have used mean fields Q — the Hartree-Fock approximation and the mean 
] field approximation for slave bosons. The 1/N correction to the latter explains the overall phase diagram rather well. 

I I The cooperative effect has been discussed on phenomenological grounds. Q Here we take a different approach, which 
'"O ' qualitatively reproduces exactly known, ground-state properties near the metal-insulator transition in one dimension. 
Ch , For M=Ag, the Ag d level is far away from the Fermi energy so that the DCNQI tt bands are quarter filled. In this 
O > sense, they are similar to (TMTTF)2X and (TMTSF)2X salts, which are also quarter-filled and show various phases. 

In the latter, the ikp anion potential has been considered to produce the umklapp process and then to cause 
the metal-insulator transition. |^ Such an extrinsic potential is absent in the (i?ii?2-DCNQI)2Ag salts, but the metal- 
insulator transition occurs, accompanied by the formation of a Akp CDW for both of the i?i=i?2= CH3 (denoted by 
DMe hereafter) and _Ri=_R2=I (denoted by DI hereafter) cases. |^ At low temperatures, (DMe-DCNQI)2Ag becomes 
' a spin-Peierls state, while (DI-DCNQI)2Ag becomes an antiferromagnet. Thus electron correlation in the tt band 
also plays an essential role to determine the ground-state phases. 

Electron correlation is essential in low dimensions. In pure one dimension, the Fermi liquid is unstable against 
any perturbation. When the perturbation is irrelevant in terms of a renormalization group, the excitation spectrum 
is gapless. Then the low-energy limit of the property is described by the Tomonaga-Luttinger liquid theory and 



o 

i> : 

0^ . 



X 



characterized by power laws in the density of states and various correlation functions. |11 Otherwise, the spectrum 
has a gap so that some of the correlation functions decay exponentially. p3[ 

For the metallic (DMe-DCNQI)2Cu, the single-particle spectrum has been observed in photoemission experiments 
and described by a power law of the electron binding energy. This suggests that the metallic phase is described 
by the Tomonaga-Luttinger liquid theory for low but not extremely low temperatures. The exponent is much larger 
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than the calculated one for the Hubbard model or the extended Hubbard model with on-site and nearest-neighbor 
repulsions only. However, it can be explained by long-rage interactions. |]l5| , |l6| The power-law holds up to the binding 
energy as large as 0.3 eV. This fact would also indicate a long-range interaction. Hi Recently, suppression of the 
interchain coherent hopping by such a large exponent in the single-particle spectrum is numerically investigated. 

In the insulator phase on the other hand, the Cu d electrons are almost localized so that good one-dimensionality 
is expected. Since the nesting is perfect and the tendency to an insulator phase is strong in one dimension when 
coupled with phonons, the metal-insulator transition must be explained at least in a purely one-dimensional model 
in order for it to occur in a quasi-one-dimensional real system. We then employ a purely one-dimensional model to 
study conditions for the observed Mott insulator phase with a gapless spin mode. 

Meanwhile, (DCNQI)2Ag salts are insulators for low temperatures. One-dimensionality is expected to be better 
than (DCNQI)2Cu salts because the metal ion does not help the electron propagation perpendicular to the most 
conducting direction as in the Cu salt. In fact, the one-dimensional band structure has been evidenced by the 
polarized reflectance spectra on single crystals. |17 Thus the Tomonaga-Luttinger hquid theory would be a good 
starting point. Qualitative aspects of the experimental findings are in fact explained by considering the effects of 
electron-2fci? phonon coupling, electron-4fci^ phonon coupling, and slight three-dimensionality on a one-dimensional 
continuum model. 

In one dimension, a renormalization-group method based on the scaling law is very useful to study the effects 
of various perturbations. The Mott transition is caused by the umklapp process. The temperature and frequency 
dependence of the conductivity for commensurate and nearly commensurate fillings has been studied in detail [ |l8| with 
the renormalization-group method combined with a memory-function approximation, [ p^ The Mott transition which 
occurs when the correlation strength is changed and that which occurs when the filling is changed have different 



physical properties, which are clarified for even and odd commensurabilities. |20| Meanwhile, the electron-phonon 
interaction produces the retarded attraction. Its effect has also been studied with the renormalization-group method. 
When it contributes to the forward scattering, the backward scattering, and the umklapp process, it enhances the 
pairing correlation, the formation of a spin gap, ||2|] and the formation of a charge gap, |^,^ respectively. 

Here, we apply the method to more complex systems, which have couplings with 2kF and Akp phonons with one or 
two bands at even or odd commensurability. We derive and numerically solve the lowest-order equations to learn which 
phase the Tomonaga-Luttinger liquid approaches as the energy scale is lowered. Although a precise description of the 
electronic states at strong-coupling fixed points is beyond the scope of the present method, it indicates whether each 
excitation spectrum has a gap or not so that the ground-state phases can be classified accordingly. Such procedure 
reproduces qualitative tendency to the strong-coupling phases for nearly-half-filled electron-phonon systems 23 
and for a two-coupled-chain system. ||2^ The present paper also follows the procedure. 

To compare with experimental results, it would be necessary to take weak three-dimensionality into account. The 
scaling law and the consequent power-law behavior of various quantities in one dimension result from the interference 
of the logarithmic singularity in the 2k f particle-hole channel with that in the particle-particle channel. The con- 
tributions from the corresponding lowest-order bubble diagrams are different in signs only. Quasi-one-dimensionality 
would be simulated by different cutoff energies for the two logarithmic singularities and the consequent imbalance 
between the two channels. In the M=Ag case, we find that the differences between the i?i=i?2=CH3 salt of the 
spin-Peierls ground state and the _Ri=i?2=I salt of the antiferromagnetic ground state are qualitatively explained by 
a difference in the cutoff energy in the particle-particle channel. 

In the M=Cu case, the two-band feature is essential. The interference between electron-phonon couplings, backward 
scatterings and umklapp processes is much more complicated in the two-band case than in the single-band case. The 
phonons with momenta near 27r/3 are responsible for the Akp scattering in the sixth-filled band and for the 2kF 
and Akp scatterings in the third-filled band, while the phonons with momenta near 7r/3 are responsible for the 2kF 
scattering in the sixth-filled band. Different combinations of the above scatterings lead to different backward scattering 
or umklapp processes, opening charge gaps in both of the sixth-filled and third-filled bands. The insulator phase can 
have a spin gap only in the third-filled band (Mott insulator) or spin gaps in both bands (Peierls insulator). The 
former phase is realized when the sixth-filled band {i.e., ir-d hybrid band P5[ ) has strong correlation and the electron- 
2kF phonon coupling is moderate. The electron-4fcF phonon coupling is necessary here, though it may not be strong. 
The latter phase is realized when the electron-2fci? phonon coupling is strong. 

To study how the behavior of resistance depends on the commensurability and model parameters, we use a memory- 
function approximation. For M=Ag, the renormalization-group method qualitatively reproduces the difference be- 
tween the DMe and DI cases. For M=Cu, the behavior of the resistance above the transition temperature is less 
sensitive to changes in correlation strengths than in the quarter-filled case. Experimentally, the transition is of first 
order due to the third-order commensurability energy, p6[ and it needs three-dimensionality in phonons, which is 
beyond the scope of the present study. 

This paper is organized as follows: Sec. || introduces bosonized models based on one-dimensi onal continuum 



models with backward scatterings, umklapp processes, and couplings with 2kF and Akp phonons. Sec. HI outlines the 
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derivation of renormalization group equations, clarifying how different electron-phonon couplings effectively produce 
different backward scattering or umklapp processes and how they open gaps in different channels. Sees. [V and 
show phase diagrams for different parameters and resistance as a function of temperature for M=Ag and M=Cu, 
respectively, to compare them with the experimental results. Sec. ^ summarizes the present work. Part of the 
results presented in this paper were reported briefly elsewhere, 



II. MODELS 



We consider a continuum model in which the noninteracting part is a Tomonaga-Luttinger liquid extended to 
include the spin (denoted by subscript cr) and charge (denoted by subscript p) degrees of freedom, 



Hn = 



v—a^p 



dx 
2^ 



(2.1) 



where fields <p,y{x) and {1/Tr)dx0^{x) are conjugate, u,^ and K^, are the velocity and the correlation exponent of the 
v channel, respectively, which are standard notations and described in detail in the previous paper. For M=Ag, 
we consider a quarter-filled band. For M=Cu, we take two bands and distinguish them by subscripts: A for the 
sixth-filled band; B for the third-filled band. Then, the noninteracting part is the sum of Haq and Hbo, in which 
the summation is performed over h'=Aa, Ap and v=Ba, Bp, respectively. It is noted that -ff(c)o (C* = A,B) [i.e.. 
Ho, Hao, or Hbo] contains the (one-electron) kinetic part and the forward scattering. For M=Cu, we consider, for 
simplicity, the regime where the Fermi velocity in the A band and that in the B band do not differ so much. Since 
we do not expect that the difference in the Fermi velocities affects the scenario to the metal-insulator transition 

vp- For later convenience, we define Ar(c)(T 



we use the averaged Fermi velocity, vp- For later convenience, we define -'^(c)ct = 2(1 — K^^^^) ~ 5(c)i/(^^f) 



and X 



{C)p = 2(1 - i^(c)p) ^ (ff(c)i - '^9(C)2)/{'^vf) for C = where g^^c)! 



^0 

and g{c)2 S'^e the backward and 
forward scattering strengths (in the C band), respectively. In terms of the single-band extended Hubbard model 
with on-site ([/) and nearest-neighbor {V) repulsions, = U/{ttvf), Xp — —{U + AV)/{ttvf) for quarter filling, 
XAa = {U + V)/iTTVF), Xap = ~{U+3V)/{ttvf) for sixth filling, and Xb^ - {U~V)/{nvF), Xsp = -{U+5V)/{ttvf) 
for third filling. 



A. Electron-electron interactions 

A backward scattering between antiparallel spins of strength ^(c)^ — 9{c)i/{'^''^f), 

H{c)a = y(c)a'i^VF^ j dx V'(c)i,s(a;)'0(c)2 -s(a;)'0(c)2,s(a;) , 



(2.2) 



where V'CCji.s ^^id 4'(c)2,s right- and left-going electrons with spin s (in the C band), is written with the phase 
field as 



H 



y(C)aVF 



dx cos 



2V20(C).(x) 



(2.3) 



where a is a cutoff parameter of the order of the inverse of the Fermi wave number. 

To study metal-insulator transition, we need to include high-order electron-electron scattering processes, which are 
to be produced in the renormalization process when electron-phonon interactions are included. In the 1 / m-filled band, 
the Fermi wave number is tt/to so that the umklapp process is written as 



Hu.Tn['^pii4'<7);yp] oc ypVpa" 



dx 



+ h.c. , 



(2.4) 



where Yp is a coupling strength such that the prefactor appears simple in the phase-field representation, a is of the 
order of the lattice spacing and set to be a, and the power should be understood as point-split hereafter. With the 
phase-field operators, it is rewritten as 



„ (ypVp)/ {'^Tra'^) J dx cos [m\/2(f>p{x)\ for even m , 

HuMM^ <P-h yp\ - I {YpVF)/{^/2^a')Sdx cos [m^/2ct> p{x)\ cos [V2(/).(x)] for odd m , 



(2.5) 
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where and hereafter less relevant terms are neglected. Thus, we have Hp = Hij^4^[(j)p;Yp] for the quarter-filled band, 
Hap = Hu,6[4>Ap]YAp] for the sixth-filled band, Hspa = Hu,^[(t>Bp,4>Ba\YBpa\ and Hsp = ^^(7,6 [0Bp ; ^Bp] (which 
is included though less relevant than Hspa) for the third- filled band. In terms of the extended Hubbard model, 
Yp = / {16TT^t^vp) — /8, etc. As m increases, the effect of Hjj^rn becomes weak on opening a gap, as explicitly 
shown in the renormalization-group equations later. Note that, for odd m, the umklapp process involves the spin 
degrees of freedom: po| when the umklapp process flows to a strong-coupling fixed point, there are gaps in both the 
charge and spin excitations. In such a case, soliton excitations carry both charge and spin in contrast to the case with 
even m. 

In addition, there are interband scattering processes for M=Cu. The most relevant backward scattering is written 

as 



HABpa = (YABpaTf VFa/V2) / dx 



h.c 



(2.6) 



where Yabp^ is the coupling strength. It is rewritten as 



HABpa = ^^"^^ J dx cos 2^2(1) Ap{x) - V2(f>Bp{x) 



cos 



V2<pBaix) 



(2.7) 



which involves the charge degrees of freedom in the A band and both the charge and the spin degrees of freedom in 
the B band. The scaling dimension is the lowest (in the noninteracting limit) among the perturbations involving the 
charge degrees of freedom. So, this is the most likely candidate which causes a metal-insulator transition. When this 
backward scattering (not umklapp process) flows to a strong-coupling fixed point, there are at least three gaps in the 
corresponding channels. Meanwhile, the most relevant umklapp process is written as 



HaBp = YABpaT^^Vpa'^ J dx 



2 r 



'^^Bl,six)^B2A^ 



h.c. 



(2.8) 



where Yabp is the coupling strength. It is rewritten as 

„ YABpVF f, 

Habp = ^^^2 ' 



2V2(l)Apix) + 2V24> Bp{x) 



(2.9) 



which involves the charge degrees of freedom in both bands but does not involve the spin degrees of freedom. 



B. Electron-phonon interactions 



As to phonons, we need to consider those with momenta near 2fci? and those with momenta near Akp. Because 
we are not interested in superconductivity here, we do not consider phonons with momenta near Okp, though the 
extension is straightforward. p^ , |2^ The wave numbers ^—2kp, Akp are /i=7r/2, tt for m=4, ^=7r/3, 27r/3 for m=6, 
and /i=27r/3, 47r/3(=— 27r/3) for m=3. The corresponding phonon fields and their conjugate momenta are denoted by 
(j)p,{x) and n^(a;), respectively. The fields (^^{x) and 11^ (x) are real, and the others are complex. The phonon parts 
of the models are written as 

Hp = j dx [Wp,{x)n^^{x) + uolcj)p,{x)^^^{x)\ (2.10) 
for the complex fields (/i ^ tt) and 

H^ = \j dx[T\l{x)+ujlcj>l{x)\ (2.11) 

for the real fields, where and lOtt are the corresponding phonon frequencies. 
Elcctron-2fci? phonon coupling is generally written as 

Hi,(c,)p. = \/7rwF'^My(c)iX! / dxip\c)2s^^)'^(C)is{x)(pp{x) -hh.c. , (2.12) 
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where y(c)i is the strength of the 2kp scattering (in the C band) by a phonon. It is rewritten as 

= '^y^^''^' f dx cos [V20(c),(x)l e-*^*«^)''(-)0^(x) + h.c. (2.13) 



For the quarter-fiUed band, for example, the couphng strength and the phonon frequencies are given by yi — 
P/y/TTVpK and = w^r = y/K/M for the Holstein couphng, [|8) J2i {(^li^i ~^ T^i ~^ 2liPi) ' ^^^^ couphng 

strength /3, spring constant K, ionic mass M, electron density rii at site i, lattice displacement qi and its con- 
jugate momentum pi, and by yi = 2ias/\/'n'VFK and \/2'^7r/2 = ^tt — 1^J~kJm for the SSH coupling, ]2g| ] 



2MPi 



ing an electron with spin s at site i, and other parameters as defined above 
Electron-4A;F phonon coupling is generally written as 



with coupling strength as, QiA+i = qi+i — qi, Q.s annihilat- 



h.c. , (2.14) 



where y[c)3 is the strength of the Akp scattering (in the C band) by a phonon. It is rewritten as 

F3,(c,). = -^^^^ / dx e-2^^(-)^(-)0,(x) + h.c. (2.15) 



27ra 

For the quarter-filled band, for example, the coupling strength is given by y^ — (3' / ^/Sir'-^vpR for the electron- 
electron-phonon coupling, '^^{U — (3' qijUi-^rii^, with on-site repulsion U, coefficient j3' of its modulation by lattice 
displacement qi, electron density Uis with spin s at site i, and by y^ = io' j y/^-K^vpK for the electron-electron-phonon 
coupling, Ei s s'^ ~ ^ qi,i+i)"^is'^i+\s' , with nearest-neighbor repulsion V , coefficient a' of its modulation by lattice 
displacement 

The phonon fields are bilinear, so that they can be integrated out completely to produce effective retarded interac- 
tions, ||22| which are lengthy and not shown here. In the ^ cx) limit, the electron-phonon interactions do nothing but 
shift the parameters, X(c)a X(c)a~ I 2/(C)i >^(C)<t Y(c)a~ I 2/(c)i X(c)p X(C)p~ I 2/(C)i P -4 | y(c)3 

Yp Yp - {yi + J/3^)/2, Yspa Yspa + {yBiyB3 + ^512/53)' ^ABpa YABpa + (J/S1J/a3 + VbiVAs), and 

Xabp ^ ^ABp - (yAsyBs + y^sJ/sa)- Below we take the Hamiltonian 

H = Hq + Ha + Hp + H^/2 + Hjr + HiT^/2 + -ffs.TT (2-16) 
for M=Ag with a quarter-filled band, and 

H = ^ (-ffcO + He a + Hcp) + Hspa + HABpa + H AB p 
C=A,B 

+HTr/3 + H2Tr/3 + Hi^a,tt/3 + HiB,2Tr/3 + -^^3,A,27r/3 + -^^3,B,-27r/3 (2-17) 

for M=Cu with coupled sixth- and third-filled bands. 



III. RENORMALIZATION EQUATIONS 

We have derived the equations following the previous study. We will give an outline for discussions later. The 
present one-dimensional quantum-mechanical system is mapped to a two-dimensional classical system where Burgers 
vectors interact with one another. The correlation function 

^j.^g»v^0(c,„(x,r)g-«\/20(c)„(o,o)^ (^C^A,B, v^a,p) (3.1) 

is perturbationally developed and successively integrated by changing the length scale little by little. In this process, 
there are two possibilities for the fate of Burgers vectors: a pair of neutral Burgers vectors are annihilated in the 
larger length scale; and a pair of non-neutral Burgers vectors are combined to produce another Burgers vector. 

2^ By comparing the correlation functions in 
successive length scales, we find relations between the effective parameters in the larger energy {i.e., smaller length) 
scale and those in the smaller energy {i.e., larger length) scale, which are described by differential equations (called 
renormalization-group equations) shown below. 



The latter is important since it causes various interference effects. |22 
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A. Quarter-filled band 



The combinations for the non- neutral Burgers vectors in the renormalization process are (jji, yl, Ya-), (2/3, ys, Yp), 
and Yp), where (a, b, c) means that any two of them are combined to produce the complex conjugate of 

the third. In order to make a comparison with the half-filled case easier, we use the notation similar to that in the 
previous study, Q =| 2/i ^3 =1 2/3 and a form factor / = (yf +y^'^)/{2Y3) which satisfies — 1 < / < 1. Below 
we mark the contributions from the particle-particle channel by Jri{l) and those from the 2kp particle-hole channel 
by Jk{1), which are useful to study a three-dimensionality effect later. Then, we finally have 

dX,{l)/dl = -J^{1) [X,{1) - Xp{l)] X,{1)I2 ~ -Ul) [X,{1) + Xpil)] X„{l)/2 - Y^{l)D^,2{l) , (3.2) 

dXp{l)/dl = Ul) [iXlil) + Xlil)] /4 - J«(0 { [iXlil) + Xlil)] /4 + AY^{1)} - Y^{l)D,,^{l) - 4Y,il)D,{l) , (3.3) 

dYp{l)/dl = Ml) [2 - SKpil)] Yp{l) - fYs{l)D^{l) , (3.4) 

dYi{l)/dl = J«(0 [2 - K,{1) - Kp{l) - X^V)] Yi{l) , (3.5) 

dY3{l)/dl = Ml) [2 - 4Kp{l) - fYp{l)] Ysil) , (3.6) 

where l=ln[EF/E], E{1) is a cutoff energy [E{Q) = Ep], Ki,{l) = (1 — Xy{l)/2)~^, is the phonon propagator 

defined by Df^{l) = [w^/£'(^)] exp[— a;^/i?(Z)]. We have used the relation Xa{l) = Yc^{l) due to the spin-rotational 
symmetry and omitted the equations for w^(/) which does not couple with the above equations in the present, 
lowest order. Initial conditions at / = are determined by treating the effective retarded interactions carefully, p^ ] 
For Xo-(O), Xp(0), and Yp{Q), only the phonon propagator is integrated from —00 to with the fixed prefactor: 
X^(0) =Xa-Yi{l- e-'^-/2/^f ), for example. 

The reason why X^ does not renormalize ^^(0) or why Yi is not renormalized at / < in this example has been 



already discussed. This correctly reproduces the antiadiabatic limit shown below ( 2.15 ). 

In real materials, slight three-dimensionality would manifest itself at very low temperatures where the electron 
transport in the perpendicular direction becomes coherent. The renormalization-group method would not be justified 
deep in the anisotropic three-dimensional regime, where the scaling law no longer holds. However, the scaling law would 
deteriorate little by little as that regime is approached from above. Then we can see at least how the renormalization 
flow is deflected and which phase it tends to approach within the present method. The scaling law in one dimension 
results from the interference of the 2kF particle-hole channel with the particle-particle channel. The corresponding 
lowest-order bubble diagrams are logarithmically divergent and have coefficients of equal magnitudes and different 
signs. Such interference disappears in higher dimensions because the Fermi surface does not consist of a finite 
number of points any more. Though the distortion of the Fermi surface removes the logarithmic singularity in the 
particle-hole diagram unless the nesting is perfect, the particle-particle channel generally becomes less important for 
repulsive interactions when only the particle-particle diagrams are summed infinitely. It is natural to regard quasi- 
one-dimensionality as causing imbalance between the two channels as the most important effect among all. [| It IS 
realized by different cutoff energies for the two logarithmic diagrams. 

It should be noted that the renormalization-group equations derived from the bosonized model and the mapping 
to a two-dimensional classical system and those derived directly from the original fermion model are equivalent. |ll| ] 
Renormalization of the velocities and detailed forms of cutoff functions if any {e.g., one due to a misfit parameter) are 
generally different, but they are not essential. Deriving the renormalization equations again in the second method, 
we can distinguish the contributions from the particle-particle channel and the contributions from the 2kp particle- 
hole channel, as already shown above. Then we can cut off either the particle-particle channel by the function Jri{l) 
(determined below) or the 2kF particle-hole channel by the function Jn{l). Assuming that the logarithm ln{Ep/E) 
is replaced by ln[£'|,/(£'^ -I- 7^£'|,)]/2 (7 = 77, k) in the perturbation expansions, we obtain 

J^il) = (d/dl) In [El/{E^il) + -f^El)] /2 = (I + 72e2')-i , (3.7) 

which satisfies Jo(0 = 1, J-y^oil) <C I for Z » — ln7 {i.e., E <^ jEp), and Jj^o{l) — > as / — > 00. 

The temperature dependence of the resistance shows how the system approaches the insulating/metallic ground 
state. It is a useful property to characterize the metal-insulator transition, and also a good indicator for the quality 
of theoretical approaches. Since we perform perturbative calculations above, we perform the perturbative expansion 
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for the conductivity following the study for commensurate and nearly commensurate fillings [Q with a memory- 
function formalism. ||l9[] The conductivity is given by (t(w) = {i2upKp/iT)[Lo + M((x')]~^, where the memory function 
M{ijj) is defined by A/(ti>) = uj{j\j)^/ {{j;j)o — (j; and is the retarded correlation function of the current 

operator, j{x) = \/2upKpdxdp{x)/'K. The lowest-order term in the perturbative expansion for M(w) is Af(w) ~ 
((F;F)° — {F;F)^) /{—u){j; j)o), where F is defined by F = [j,H] and stands for the retarded correlation 

function of the operator F at frequency to computed in the absence of perturbations. Later we consider a temperature 
range below the phonon frequencies, where electron-phonon interactions are integrated out. The calculation of {F; _F)2, 
is straightforward. Taking the oj ^ limit, we finally have, for the resistance R{T) = l/a{ui = 0, T) at finite 
temperature T, 

cx Y^{T)TB^ [4ifp(r), 1 - 8Kp{T)] cos^ [4^i^p(T)] , (3.8) 

where B{x) is the beta function and Yp{T) denotes Yp[l = ln(£'i?/T)], etc. If we neglect the temperature dependence of 
Kp in the equation for Yp (and for k = 0), which would be valid only in the weak-coupling limit, R{T) cx y^T^^^p^^ . 



B. Coupled sixth- and third-filled bands 

The combinations for the non- neutral Burgers vectors in the renormalization process are now (y^a, Ubi^ Yabpo-), 

ySZ, YaBp), {VBI, J/B3) ^Bpa), O^ABpa, ^ABp, ^Bpa), O^ABpa, ^ABpa, Ybct), O^Bpa, ^Bpa, Xbct), {YBpa, YBpa, 

Ybp), iVAi, y^i, YAa-)i and (ubi, Vbii ^Bct)- When (a, b, c) is possible, (a*, 6*, c*) is also possible but it is not listed 
above because it is trivial. Note that y's are complex parameters, while y's are real parameters. Finally, we have 



dXAa{l)/dl = ~~Xl^{l)- I yAi{l) P D^,^{1) , (3.9) 

dXAp{l)/dl = -YispAl) - YlBp{l) - 9Yip{l)- I yAiil) P D^/,{1) - 4 | y^3(0 P D,^/,{1) , (3.10) 

dXBa{l)ldl = ~Xl,{l) - \yIspAI) - \ybpAI)- I VBiil) ? D2./^{1) , (3.11) 

dXBp{l)/dl = ~\yIsp,{1) - Ylsp{l) ~ ^rj,,(/) - 9rJ,(/) - [| yBi{l) I' +4 I ys3(0 ?] D^^Ml) , (3.12) 

dyAi{l)/dl = 1^1 - ^[KAa{l) + KapH) + XAAl)]j VAiil) , (3.13) 

dyAz{l)/dl = [1 - 2KAp{l)] yAsil) + ^i^ABpa(/)2/Bi(0 - lYABp{l)y*B3{l) : (3-14) 

dyBl{l)/dl = {l - ^ + KBp{l) + XBail)]^ 2/51(0 + \YABpa{l)yA3{l) + ^i^Bp. (0J/B3 (0 , (3-15) 

dyB3{l)ldl = [1 - 2i^Bp(0] ys3(0 - \YABp{l)y*Az{l) + \YBp.{l)y*BS) , (3.16) 

dYABpa{l)ldl = {2 - ^ [4^Ap(0 + i^S.(0 + KBp{l) + Xi3.(0]| i^ABp.(0 

- \YABpmBpAl) + [yBl(0j/l3(0 + y*Bl{l)yAm i?2V3(0 , (3.17) 

dYABp{l)/dl = {2 - 2 [ifAp(0 i^Bp(0]} ^ABp(0 

- \YABpamBpAl) - [yA3(0yi33(0 + 2/A3(0yS3(0] ^2^/3(0 > (3.18) 



7 



dYBpa{l)/dl =\j-\ [KBa{l) + 9Kbp{1) + XBa{l) + Ybp{1)]^ YBpa{l) 

- IyabpAI)Yabp{1) + [yBi{i)yB3{i) + y*Bi{i)y*B3{i)] D^^,^{i) , (3.19) 

dYAp{l)/dl = [2 - 18KAp{l)] Yap{1) , (3.20) 

dYBp{l)/dl = [2 - ISKbpQ)] YbpQ) - \y^,AI) > (3-21) 

where the cutoff energy and the phonon propagators are defined as before, Xca{l) = Yca-{1) {C = A,B) has been 
used. Initial conditions at / = are also determined as before. 

For the present system, the resistance R{T) in a temperature range below the phonon frequencies is given by 

RiT) (X 5Y^BpAT)T cos^ [HKap{T) + Kbp{T) + KBa{T)}TT/A] 

X B2 [{AKap{T) + Kbp{T) + KBa{T)} /4, 1 - {4KAp{T) + Kbp{T) + KBa{T)} /2] 

+ 8YiBp{T)T cos^ [{KAp{T) + KBp{T)) tt] [ KAp{T) + i^s,(T), 1 - 2 {KAp{T) + i^Bp(T)}] 

+ 9y|^,(T)Tcos2 [{QKbp{T) + Kb,{T)} -k/A] [{QKbp{T) + Kb^T)} /4, 1 - {9i^s,(T) + KbAT)} /2] 

+ 36Y:|^(T)Tcos2 [9KAp{T)n] B^ [9KApiT), 1 - 18Kap{T)] 

+ 36rJ^(r)Tcos2 [9KBpiT)Tr] B^ [9KBpiT), 1 - 18Kbp(T)] . (3.22) 

In the above formula, any perturbation which contains (j)Ap or (j)Bp contributes to the resistance. Among them, 
Yabpct has the lowest scaling dimension and dominates the resistance. In the weak-coupling limit, we have R{T) oc 
Y^gp^T^^'^f'~^^^''~^^^>'~^ . The actual temperature dependence is modified from this simple power law. 



IV. RESULTS FOR QUARTER-FILLED BAND 

The fixed points of the spin and charge correlation exponents, K* and K*, determine the asymptotic correlation 
functions. Assuming a gap in the v spectrum for K* = 0, we classify the groimd-statc phases into a gapless state 
in which all perturbations are irrelevant so that the low-energy limit is regarded as a Tomonaga-Luttinger liquid 
(denoted by TL in the phase diagrams below), a state with only a spin gap which we call a spin-gap state (denoted by 
SG), a state with only a charge gap which we call a Mott insulator (denoted by MI), and a state with both spin and 
charge gaps which we call a spin-Peierls state (denoted by SP). When we consider three-dimensionality, the gapless 
(thus metallic) state does not correspond to a Tomonaga-Luttinger liquid any more so that it is simply denoted by 
M. The Tomonaga-Luttinger liquid has a dominant 2kF SDW correlation exp[—{K* + l)lnr] or a dominant Akp 
CDW correlation ~ exp[— 4iC* Inr]. The spin-gap state has a dominant 2kF CDW correlation ~ exp[—K* Inr]. The 
Mott insulator has a dominant Akp CDW correlation. 

Usually a spin-Peierls state is described by a localized spin system since the charge gap is much larger than the 
spin gap. Here we call even a state with comparable magnitudes of gaps a spin-Peierls state because it is smoothly 
connected with the so-called spin-Pcicrls state. We assume that phonons arc also one- dimensional for simplicity so 
that we do not consider the possibility for a static lattice distortion. In any case, a coupling with 2kF phonons is a 
necessary condition here for a spin gap and thus for a spin-Peierls state because we consider repulsive electron-electron 
interactions only. With a coupling with three-dimensional phonons, the spin-Peierls state here would be accompanied 
by a static lattice distortion. 

Before we use parameters for (DCNQI)2Ag, wc first show how the clectron-4fci;- phonon coupling and the umklapp 
process interfere with each other and next show how they are affected by the electron-2fcF phonon coupling, large 
phonon frequency and nearest-neighbor repulsion. For the bare phonon dispersion, we take a;^/2 = ^^tt/VS. It is noted 
that the quantity / appears with Yp(l) as a prefactor in the rcnormalization equations so that we fix / = — 1 and 
allow Yp to be negative. We do not cut off the logarithmic singularities, k = r] = 0, unless explicitly mentioned. 



A. Coupling with 4kF phonons and umklapp process 

It is noted again that the perturbative expansion gives Yp = /8 for the extended Hubbard model. The umklapp 
process is of third order with respect to electron-electron interaction, so that it is generally weak. In this and the 
next subsections, we vary the strength of the umklapp process to see the interference effect. 
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To have a char ge g ap, if* must vanish so that Xp{l) must approach — oo in the I — > oo hmit. From the renormal- 
ization equation (|3.3| ), it requires Yp{l) not to vanish. Then from {^A), Kp{l) must be smaller than 1/4 [Xp{l) must 
be smaller than —6] for large I. For the pure Hubbard model, it is known that it cannot be smaller than 1/2. ]3C| ] 
Even for the extended Hubbard model with nearest-neighbor repulsion, it is at least 3/16. jl^] Such a small Kp{l) is 
achieved generally for long-range repulsion. pl| ] The photoemission experiment [ p4[ may indeed suggest a long-range 
repulsion for (DMe-DCNQI)2Cu. In the (DCNQI)2Ag salts, the Akp CDW is observed §] and the lattice would be 
modulated. So, it is reasonable to regard electron-phonon interaction as cooperating with the umklapp process to 
have such a small value of Kp{l). This is the case when the umklapp process itself is very weak because it occurs as 
a high-order process. 

In the absence of the electron-phonon coupling, the critical value for a finite charge gap is large (Fig. |l|) and beyond 
the scope of the perturbative regime. However, with the electron-4A:i? phonon coupling, it becomes smaller. For 
Yp > corresponding to / — —1, e.g., for lattice modulation of the nearest-neighbor repulsion, the umklapp process 
is constructively interfered with the electron-4fcir phonon coupling, as is expected from the renormalization-group 
equations. This is the reason for the negative slope of the critical value of at Yp = in the figure. It is noted that, 
even for Yp < corresponding to / = -1-1, e.g., for lattice modulation of the on-site repulsion, the critical value of Y3 
does not become so large. This is because Yp{l) changes the sign after it decreases due to the destructive interference 
at the initial high-energy scale. Thus, the interference is constructive at low-energy scales irrespectively of the sign 
of /. This situation is in contrast to the half-filled and nearly-half-filled cases, where the opening of a spin gap is 
sensitive to the form factor of the electron-2A:p phonon couphng. |22 2^ 



B. Coupling with 2kp phonons 



Th e electron-2fci? phonon coupling is expected to decrease Xp{l), enhancing the tendency for a charge gap, from 
( |3.3[ ). However, its tendency is very weak as clearly seen from a comparison of Fig. |l| with Fig B an d directly from 
Fig. ^. There is no interference of the electron-2fcF phonon coupling with the umklapp process in ( p.4[) . This indicates 
again how the interference effect is important to cause the metal-insulator transition. 

The role of the electron-2A:i;' phonon coupling is mainly to enhance the tendency for a spin gap. For a sufficiently 
large coupling, the system has a finite spin gap (Fig. |^) as expected from (3.2). Then, the 2kp CDW correlation 
becomes dominant. Note the phase boundary between TL and SG and that between MI and SP depends considerably 
on the electron-4A;F phonon coupling. This is due to the fact that the clcctron-2fcF phonon coupling is affected by 
both of the spin and charge correlation exponents in (3.5). Meanwhile, the electron-4/ci? phonon coupling is affected 
by only the charge correlation exponent in ( [3.6| ). E ven in the weak-coupling limit [Ka{l) ~ K n(l) — 1], Yi{l) is largely 
affected by the correlation exponents as shown in (3^), although Yz{l) always decreases in ( |3.6| ). Thus, the opening 
of a spin gap is rather sensitive to the electron-4fci? phonon coupling. 



C. Phonon frequency and nearest-neighbor repulsion 

Now we have a clear idea about how the electron-phonon couplings interfere with the electron-electron interaction 
and how they determine the ground-state phase. In this subsection, we show how the phase diagram depends on 
the phonon frequency and the relative magnitude of | Xp \ to Xa, i.e., the nearest-neighbor repulsion in terms of 
the extended Hubbard model. Note that -Xp/X^ = 1 + W/U in this model. For V/U =1/4 and 1/2, the relative 
magnitude of | Xp \ becomes twice and three times as large as that for the Hubbard model, respectively, so that the 
results can be largely affected by the nearest-neighbor repulsion strength. 

The phonon frequency relative to the Fermi energy controls the energy scale of the interference. In the antiadiabatic 
limit, the electron-phonon coupling simply shifts the strengths of the electron-electron scattering parameters. For finite 
phonon frequency below the Fermi energy, the effective interaction is retarded and the interference effect becomes 
l arg est at an energy scale comparable to the phonon frequency as in the half-filled and nearly-half-filled systems. 
|P2| , P3| As the phonon frequency increases, the critical coupling strength of Y3 for opening a charge gap decreases, 
while that of Yi for opening a spin gap increases [Figs, ^(a) and (b)]. Then the Mott insulator phase occupies a wider 
area in the phase diagram. 

As the magnitude of | Xn \ increases, the initial value of Kp decreases. Then the critical value of Y3 for opening a 
charge gap decreases [Fig. Q(c)]. The slope is steeper for smaller phonon frequencies. Meanwhile, the Yi coordinate 
of the cross section of the phase boundaries among TL, MI, SG, and SP does not depend so much on | Xp \. 
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D. Phase diagram for (DCNQI)2Ag 



According to the local density functional theory the band width for (DMe-DCNQI)2Ag is about 0.9eV. Phonons 
due to the dimerization are observed at about O.OSeV in the infrared spectrum |3|]. We use •\/2a;^/2 = o;^ = OAEp- 
The condition for the Mott insulator phase is easily achieved because the phonon frequency relative to the Fermi 
energy is not small and because the nearest-neighbor repulsion relative to the on-site repulsion is not small either. 
p3[ We show the data again for completeness (Fig. ||). 

The author does not know the electron-electron and electron-phonon interaction strengths, but the qualitative aspect 
of the present results would not change for different electron-electron interaction strengths. Since (DMe-DCNQI)2Ag 
becomes a spin-Peierls state at low temperatures (in the sense that both of the spin and charge gaps open) , we expect 
that its electron-phonon coupling strengths are in the SP phase of the figure. Note that all of the four possible phases 
appear in the phase diagram in the purely one-dimensional case. This is because the Tomonaga-Luttinger liquid phase 
is stable for weak electron-phonon couplings. We consider a cutoff in the logarithmic singularity below, where the 
Tomonaga-Luttinger liquid is not realized. Then some phases disappear from the phase diagram. 



E. Effects of a cutoff in the logarithmic singularity 



We consider three-dimensionality by cutting off the logarithmic singularity in either the particle-particle channel or 
the 2kp particle-hole channel as explai ned in the previous section. The most important effect appears in the charge 
degrees of freedom. It is obvious from (|]^) that a cutoff in the particle-particle channel, ij 0, leads K* to 0, while 
a cutoff in the 2kp particle-hole channel, k 7^ 0, leads K* to infinity. Therefore, a charge gap opens in the former 
case, while it does not open in the latter case. Then, possible phases are a Mott insulator and a spin-Peierls state in 
the former case (Fig ||), and a gapless metallic state and a spin-gap state in the latter case (Fig 

If the nesting property is maintained in the quasi-one-dimensional case, the particle-particle channel becomes less 
important at low energies. Then, it is reasonable that the nesting causes a finite charge gap or makes the ch arge gap 
larger. Note that the phase boundary for a finite spin gap is shifted to the right in Fig. |^, as expected from (3.2). As 
the cutoff increases, a transition occurs from the spin-Peierls state to the Mott insulator at zero temperature. This 
is consistent with the experimental results. (DMe-DCNQI)2Ag has strong anisotropy so that it is regarded as a good 
one-dimensional material and it becomes an insulator below about 120 to 150K. Meanwhile, (DI-DCNQI)2Ag has a 
considerable conductivity in the transverse direction and it is an insulator already at room temperature. From the 
activation plot, the charge gap is estimated to be 490K ||l^. (DMe-DCNQI)2Ag becomes a spin-Peierls state at about 
80K so that it has a finite spin gap at zero temperature. Meanwhile, (DI-DCNQI)2Ag becomes antiferromagnetic 
below 5.5K so that the spin excitation spectrum is gapless. The Mott insulator in the present study is expected to 
become an antiferromagnetic [or spin-density-wave (SDW)] state when weak three-dimensionality is taken into account 
because the repulsive interaction would produce an effective antiferromagnetic coupling in the transverse direction. 
The above result is reminiscent of the earlier work for localized spins at half filling which showed the instability of 
the spin-Peierls state in quasi-one dimension against a SDW. 

If the nesting property is lost in the quasi-one-dimensional case, the possibility for a finite charge gap would be 
very strong electron correlation, which is beyond the scope of the perturbative renormalization-group approach. This 
possibility may not be excluded, but we can say at least whether the temperature dependence of the resistance 
within the present approach is qualitatively consistent with the experimental data or not. It will be done in the next 
subsection. The weakened tendency for a spin gap is common with both cutoffs. Note that the tendency for a spin 
gap is substantially suppressed here (Fig. k ^ 0) compared with the other case (Fig. |[ ?y ^ 0). This is reasonable 
in that, in the three-dimensional case free from nesting (k ^ 0), a spin gap opens when the phonon- mediated effective 
attraction overcomes the repulsion, while such condition is not necessary for the spin gap in a spin-Peierls state. 



F. Temperature dependence of the resistance 

Our previous work [ p7| has shown results which are not yet converged so that we show the convergent results here. 
Here, temperatures are in the unit of Ep ~ 2400K. When the particle-particle channel is cut off, the resistance increases 
and the metal-insulator transition temperature increases (Fig. |^). With a slight three-dimensional component, the 
resistance increases rapidly below the transition temperature. This is obtained by renormalizing both Yp and Kp. If 
only Yp is renormalized, the resistance is not so steep below the transition temperature. Including this, the overall 
and qualitative behavior of the resistance is indeed consistent with the experimental one for (DMe-DCNQI)2Ag and 
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(DI-DCNQI)2Ag. Thus, we consider that a cutoff in the particle-particle channel imitates the effect of three- 
dimensionality for (DCNQI)2Ag salts. When the 2kp particle-hole channel is cut off on the other hand, the resistance 
decreases and the system finally becomes metallic at low temperatures (Fig. This is obviously contradictory with 
the experimental data, as it is expected from the renormalization-group equations. 

The real difference between (DMe-DCNQI)2Ag and (DI-DCNQI)2Ag would not be limited to the difference in the 
three-dimensionality. In fact, (DI-DCNQI)2Ag has smaller bandwidth j^] so that the relative strength of electron- 
electron interaction would be larger than (DMe-DCNQI)2Ag. We calculated the resistance with increasing coupling 
strengths and found that the behavior is similar to the case with increasing cutoffs in the particle-particle channel 
(not shown). However, the spin degrees of freedom is almost unaffected. Therefore, the three-dimensionality is the 
key to understand their differences, although the relative coupling strength is also different and contributes to the 
difference in the resistance. 



V. RESULTS FOR COUPLED SIXTH- AND THIRD-FILLED BANDS 



The ground-state phases can be classified according to the channel(s) whose excitation spectrum has a gap. When 
any perturbation is irrelevant and vanishes in the low-energy limit, the limit is described as a Tomonaga-Luttinger 
liquid without any gap. When all of the four ATj^'s diverge and both bands have spin and charge gaps, the system 
would be basically described as a band insulator (if the spin and charge gaps are comparable in the magnitude). It 
would be realized if a CDW of period 6 is formed and the Brillouin zone is folded at ±7r/6 so that both bands have 
a gap at the Fermi points. Then, it may be called a Peierls insulator because it can occur by the Peierls mechanism 
only, i.e., without the help of electron correlation. The experimentally observed, insulator phase have a CDW of 
period 3. There are charge gaps in both bands and a spin gap in the third-filled band only. In fact, it has a long- 
ranged antiferromagnetic order at low temperature. Since it can not be described by a band picture, we call it a 
Mott insulator. In the phase diagrams below, each phase is denoted by the number of gap(s). The three important 
phases above, the Tomonaga-Luttinger liquid, the Mott insulator, and the Peierls insulator, are denoted by 0, 3, and 
4, respectively. 

There are also other phases which do not generally occupy a wide area in the phase diagrams below. They can 
be artifacts of the present, lowest-order renormalization-group approach if the occupied area is very narrow. The 
mechanism of each phase is discussed below. The phase 1 has a spin gap either in the B band (if the correlation in 
the A band is comparable to or stronger than that in the B band) or in the A band (if the correlation in the B band 
is much stronger). The former would correspond to a freely moving (unpinned) CDW of period 3. The phase 2 has 
spin gaps in both bands (if yi ^ 0) or charge gaps in both bands (if yi =0). The former would correspond to a freely 
moving, CDW of period 6. 



A. Interference between different interactions 



Before looking at the numerical results, we should consider what is suggested by the equations. When the fixed point 
of the correlation exponent K* is zero, th ere is a gap i n the v channel. For that, X^(Z) must diverge to — oo as I goes to 
oo in the corresponding equation among (^^) - ( |3.12 ). Since the phonon propagators exponentially decrease and thus 
the electron-phonon interactions (j/'s) are finally integrated out, the behavior of the electronic perturbations, YABpa{l), 
Yabp{1), Ybpit{1), Yap{1), and Ybp{1), is of particular interest. If some of them diverge, they make the corresponding 
X^{1) diverge to — oo in Eqs. ( |3.9D - ( |3.12 ). The behavior of each perturbation Y{1) largely depends upon the factor 
multipHed by Y{1) on the right-hand side of dY{l)/dl [for example, 2 - i [AKapH) + Ksail) + Kspil) + Xsail)] for 
YABpail)]- Note that, before the equations for Yao-(0 ^^'^ ^Bo-(0 ('^ot shown) are linearized to retain the spin-rotational 
symmetry, the factors were 2 — 2KAcr{l) for Y^ail) and 2 — 2Kb„{1) for Ysail)- If we neglect the interference effect, i.e., 
if we do not consider the renormalization process where a pair of non-neutral Burgers vectors are combined to produce 
another Burgers vector, these factors are given by (/ dependence is implicit hereafter) 2 — j{AKap + + Kbp) for 
YABpa, 2 - 2{Kap + Kbp) for Yabp, 2 - \{KBa + ^Kbp) for Yepa. 2 - ISKap for Yap, 2 - ISKbp for Ygp, 2 - 2KAa 
for YAa, and 2 — 2KBa for YBa- They are 2 — "scaling dimensions" in the field-theoretical terminology. As the scaling 
dimension becomes lower, the perturbation generally becomes more relevant. In the weak-coupling limit (X^ —> 
and Ki, — 1), they are —1 for YABpa, —2 for Yabp, —3 for Yspa, —16 for Yap and Yep, and for Yact and Ybct- The 
corresponding factors for the electron-phonon interactions are for yAi and ysi and —1 for y^s and yss- Namely, 
Yao-, Ysa, and yBi are marginal and the others are irrelevant in this limit. Then, a charge gap would not open 
in either band. 
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Of course, the relevance depends upon the correlation strengths. As ii'y's become smaller, the perturbations 
become more relevant. More importantly, the interference effects largely affect the relevance of each perturbation in 
the renormalization process. The strengths of Yabpct, Yabp, Yspa, Yap, and Yep are initially very small (set to be zero 
in the numerical calculations in the next section) because they correspond to high-order electron-electron scattering 
processes. However, the electron-27r/3 phonon interactions, yA3, Ubi, and ysa, interfere with one another to produce 
YABpa, Yabp, Yspa, and the latter three then interfere with one another as well as with the electron-27r/3 phonon 
interactions. Once YABpa becomes relevant and diverges, Xap, Xsa, and Xsp diverges to — oo and gaps open in the 
corresponding channels. If Yabp becomes relevant, gaps open in the channels Ap and Bp. If Yspcr becomes relevant, 
gaps open in the channels Ba and Bp. Since the scaling dimension of YABpa is the lowest among the three, the first 
situation would always occur when the second or the third situation is realized. Without the electron-4A:ir phonon 
interactions, uas and i/bs, on the other hand, YABpa, Yabp, Yspa, and Ysp remain zero (when they are initially set 
to be zero), so that a charge gap does not open in either band. In such a case, the situation for yBi and XBa is the 
same as that for yAi and XAa discussed below. In contrast to YABpa, Yabp, YBpa, and Fsp, the perturbation Yap 
does not interfere with any other perturbation. Its scaling dimension is much higher than the others so that it would 
not play an important role. 

It is noted that, if YABpa or YBpa is finite (though it may be very small in the real materials), it and yBi would 
interfere and produce yA3 or yBs, respectively. This is also a possibility for opening charge gaps if an infinitesimal 
(or very small) value of yAs or yBs is necessary in the numerical results with a vanishing (initial value of) YABpa 
or YBpa in the next section. In the present two-band system, it is essential that phonons with momenta near 27r/3 
contribute to yAi, yA3, and yB3 processes. In the single-band quarter-filled case, the coupling with 2kF phonons 
interferes with the backward scattering, and the coupling with Akp phonons interferes with the umklapp process, but 
the two couplings do not interfere with each other(Sec. Ill A) in contrast to the present case. 

The electron- 7r/3 phonon interaction yAi does not interfere with another electron-phonon interaction. It only 
interferes with XAa- After yAi is integrated out, i.e., \ yAi P ^7r/3 vanishes a.t I > Iq for some Iq, XAa follows 
XAail) X^^Ho)]-'^. The fixed point is either X\^ = [K\^ = 1] or X\^ = -oo [K\^ = 0], depending 

upon the sign of ^^^^(^o)- In other words, when the electron-7r/3 phonon interaction yAi is strong enough for ^^^^(^o) 
to be negative, a gap opens in the channel Aa. Otherwise, the excitation spectrum of this channel is gapless. If 
the electron-4/ci? phonon interactions, yA3 and yB3, are absent, the behavior of yBi and XBa is exactly the same 
as that of yAi and XAa mentioned above. Without the electron-2fc^ phonon interactions, yAi and yBi, YABpa and 
YBpa remain zero (when they are initially set to be zero), so that a spin gap does not open in either band for the 
repulsive case, where the initial conditions for XAa and XBa are positive. In short, electron-2fci? phonon coupling is 
necessary for a spin gap, while electron-4fci? phonon coupling is necessary for a charge gap unless a sufficiently strong 
long-range interaction makes the scaling dimension of a high-order electron-electron scattering process lower than 2. 
This statement holds also for the single-band quarter-filled case. (Sec. iV) 



B. Phase diagrams 

In the phase diagrams below, the charge gaps in the phases 3 and 4 are brought about by YABpa, while the spin 
gap in the A band in the phases 2 and 4 are brought about by yAi- The spin gap in the B band is caused by either 
YABpa or yBi. We initially set yAi = 2/bi = 2/i and yA3 = ys3 = 2/3 and vary yi and j/s for several sets of {XAa, Xap, 
XBa, Xbp), where Xap and Xbp are negative. The strengths of Yabpo-, Yabp, YBpa, Yap, and Yep are initially {i.e., 
at / = — oo) set to be zero in the numerical calculations since they are of high order with respect to electron-electron 
scatterings. For the bare phonon dispersion, we take 2w7r/3 = /3/ = 0.4£'i? as in the previous section unless 
explicitly mentioned. 



When electron correlation is weak in both bands (Fig. 10), the phase diagram mainly consists of the phase 
(Tomonaga-Luttinger liquid) and the phase 4 (Peierls insulator), as expected, unless only the electron-4fci? phonon 
coupling 2/3 is strong. As the electron-2fci? phonon coupling yi increases, a spin gap opens in the B band either mainly 
by yBi (phases — > 1) or by constructive interference of j/si with yA3 and yB3 (phases — s- 3). The phase 1 may be 
changed into the phase 3. As yi further increases, a spin gap opens also in the A band (phases 1^2, phases 3^4). 
It should be noted that the correlation strength {XAa, Xap), the electron-2fci? phonon couphng strength {yAi), and 
the electron-4fci;' phonon coupling strength {yA3) in the A band are the same as the corresponding coupling strengths 
in the B band in this figure, but a spin gap opens first in the B band. This is due to the constructive interference of 
ysi with yA3 and yB3- In fact, the critical coupling strengths for spin gaps are the same for 2/3 = and different for 
2/3 7^ 0. As j/3 increases, the difference becomes large because the critical coupling strength for a spin gap in the B 
band becomes small. 



When electron correlation is strong in the A band (Fig. 11), a stronger electron-2A:i? phonon coupling yi is necessary 
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for a spin gap in the A band, as expected. If the correlation remains weak in the B band, the relation between the 
critical coupling strength for a spin gap in the B band and that for charge gaps is very similar to the weakly correlated 
case. Namely, as yi increases, a spin gap opens in the B band either mainly by ysi (phases — > 1) or by constructive 
interference of ysi with yAs and yss (phases — *■ 3). Then, the phase 1 is changed into the phase 3. Finally, when yi 
overcomes the strong correlation in the A band, a spin gap opens in the A band (phases 3 — > 4). As a consequence, 
the phase diagram mainly consists of the phases (Tomonaga-Luttinger liquid), 3 (Mott insulator with a CDW of 
period 3), and 4 (Peierls insulator with a CDW of period 6). It should be noted again that the electron-4A;F phonon 
coupling is necessary for charge gaps. On the j/i axis (2/3 = 0), the possible phases are 0, 1, and 2. Therefore, the 
condition for the phase 3 is strong correlation in the A band, moderate electron-phonon couplings (not too strong to 
overcome the electron correlation in the A band), and finite electron-4A:i? phonon coupling though it may be small. 



On the other hand, when electron correlation is strong in the B band (Fig. 12), a stronger electron-2fci? phonon 
coupling yi is necessary for a spin gap in the B band, unless ys is strong. Meanwhile, the critical coupling strength 
for a spin gap in the A band remains small if the correlation remains weak in the A band. This situation does not 
correspond to (DCNQI)2Cu salts because the correlation is expected to be stronger in the A band. If 2/3 is small, a 
spin gap in the A band opens first and then the other three gaps open as yi increases. If y3 is large on the other 
hand, the three gaps except the spin gap in the A band open first. Note the phase 1 here is different from that in the 
other figures in the sense that the A band has a spin gap here. 

When electron correlation is strong in both bands and their strengths are the same (Fig. 13), the phase 1 disappears 
and the phase 2 exists only if j/3 = (spin gaps) or if yi — (charge gaps). The critical coupling strength for a spin 
gap in the B band is the same as that in the A band for 2/3 = and becomes smaller than the latter once 2/3 7^ as 
in Fig. |o|. The two curves for these critical coupling strengths are shifted to the right for small 2/3, compared with 
Fig. [^. 

So far we took Xao- =\ \ and XBa =| Xsp |, assuming only the on-site repulsion in electron-electron inter- 
actions. When we consider nearest-neighbor repulsion, for example, such a relation no longer holds. Recall that, in 
terms of the single-band extended Hubbard model, Xau — [U + V)/{itvf)^ Xap — ~{U + 3V)/{nvF) for sixth fiUing, 
and Xbo- = [U — V) / {-Kvp), Xsp = —{U -\- 5V) / {ttvf) for third filling. Even for the same coupling strengths in terms 
of the extended Hubbard model, the scattering strengths depend on the filling factor. Neglecting the difference in the 
Fermi velocities for simplicity and taking the X values of the single-band extended Hubbard model for sixth and third 
fillings, we study the effect of the nearest-neighbor repulsion with U /{nvp) = 0.4 and V/iTrvp) = 0.1 (Fig. [I4[ ). The 
resultant phase diagram is similar to Fig. ^because Xact > Xsa- When compared with the result for U/{Trvp) — 0.4 
and V/^TTVp) = (Fig. the nearest- neighbor repulsion increases the critical coupling strength for a spin gap in 
the A band and decreases that for a spin gap in the B band, as expected from the fact that the nearest-neighbor 
repulsion increases Xao- and decreases Xbo-- 

Experimentally, pressure induces the metal-insulator transition. The transition is of first order due to the third- 
order commensurability energy | p6| and it needs three-dimensionality in phonons, which is beyond the scope of 
the present study. As temperature decreases, the resistance abruptly increases at the transition under pressure for 
(i?ii?2-DCNQI)2Cu with i?i=i?2=CH3 or I, but above the transition temperature the resistance behaves as at ambient 
pressure where no transition occurs. The pressure would change the band width at least, but it hardly affects the 
behavior of the resistance above the transition temperature. We decrease the band width by multiplying X's, \ y p's, 
and uj's (all of which scale as the inverse of the band width) by a common factor (1.1 to 1.4) and calculated the 
temperature dependence of the resistance in the case of Fig. |ll| with 2/1 and 2/3 near the boundary between the 
phases and 3 (Fig. ^5|). The overall behavior of the resistance is insensitive to changes in the parameters even near 
the metal-insulator transition, which is qualitatively consistent with the experimental results. This is in contrast to 
the quarter-filled case appropriate for (DCNQI)2Ag salts, where the correlation strength affects the behavior of the 
resistance rather sensitively. (Sec. IV) 



VI. SUMMARY 



Metal-insulator transitions and electronic phases in (DCNQI)2M (M=Ag, Cu) salts are studied with the 
renormalization-group method for the one-dimcnsional continuum models with backward scatterings, umklapp pro- 
cesses and couplings with 2kp and Akp phonons. These salts are in contrast to the quarter-filled (TMTTF)2X and 
(TMTSF)2X salts, where the extrinsic Akp anion potential produces the umklapp process. For the present salts with 
M=Ag, such a potential is absent but the electron-4fci? phonon coupling interferes constructively with the umklapp 
process, thereby causing a metal- insulator transition. The Akp CDW is therefore a product of the cooperation of 
the electron-electron and electron-phonon interactions. It can be viewed as a Mott insulator in the sense that it has 
gapless spin excitations and it cannot be described by a band picture. 
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Experimentally, the physical properties of (i?ii?2-DCNQI)2Ag depend upon 7?ii?2- For _Ri=i?2=CH3, one- 
dimensionality is rather good. It becomes first an insulator at about 120K and then a spin-Peierls state at 80K, 
opening a spin gap. Meanwhile, for Ri=R2=l, electron transfer in the perpendicular direction is not negligible and 
electron correlation is expected to be stronger due to the narrower bandwidth. It is already an insulator at room tem- 
perature and the spin excitation spectrum remains gapless to zero temperature. In fact, it becomes antiferromagnetic 
below 5.5K. Such qualitative difference is explained by the present approach if the three-dimensionality is taken into 
account by cutting off the logarithmic singularity in the particle-particle channel. The cutoff suppresses the opening of 
a spin gap and enhances the charge gap. The stronger correlation for Ri=R2=l would contribute to further enhancing 
the charge gap. The temperature dependence of the resistance is calculated with a memory-function approximation 
and its behavior is consistent with the experimentally observed one. 

For M=Cu with coupled sixth- and third-filled bands, the Mott transition is accompanied by the formation of a 
CDW. In real materials, three-dimensionality would not be neglected. Reasons why we take the one-dimensional 
model arc: a Tomonaga-Luttingcr behavior is observed in the metallic phase in photoemission experiments; and the 
metal-insulator transition should be explained in one dimension also. 

In order for gaps to open in the charge excitations in both bands and in the spin excitation in the third-filled 
band, the high-order backward scattering YABpa must be relevant. Its scaling dimension is higher than 2 in the 
noninteracting limit. So, it may look as an irrelevant perturbation. It is not the case if the interference of the 
electron-2A;F phonon coupling ysi with the electron-4A:F phonon couplings, yAs and yB3, is taken into account. 
These electron-phonon couplings interfere with one another and produce effective high-order backward and umklapp 
scatterings. The interference of the electron-phonon couplings and the high-order scatterings is constructive and can 
make Yabpct a relevant perturbation. 

In order for a gap not to open in the spin excitation in the sixth-filled band, the electron correlation must be strong 
enough. Therefore, the condition for the experimentally observed, Mott insulator phase with a CDW of period 3, 
i.e., with a gapless spin mode, is strong correlation in the sixth-filled, n-d hybrid band, moderate electron-phonon 
couplings which are not too strong to overcome the strong electron correlation in the sixth-filled band, and finite 
electron-4fcF phonon coupling which may be small. The temperature dependence of the resistance is found to be 
insensitive to changes in the parameters even near the metal-insulator transition, which is again consistent with the 
experimentally observed one. If electron-phonon couplings were too strong, the insulator phase would be accompanied 
by a CDW of period 6 and have gaps in all channels. 

Showing these results, we have demonstrated that the renormalization-group method reproduces the qualitative 
aspects of the ground-state phases and the behavior of the resistance in these quasi-one-dimensional organic materials 
very well. 
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FIG. 1. Phase diagram for = 0.4, Xp = -0.8, = O.lEp, and Yi = 0. 



FIG. 2. Phase diagram for X^ = 0.4, Xp = -0.8, uj^ = O.lEp, and Yi = 0.2. 



FIG. 3. Phase diagram for X^ = 0.4, Xp = -0.8, Yp = 0.375X^, and = O.IEf- 

FIG. 4. Yi and Yi coordinates of the cross section of the phase boundaries among TL, MI, SG, and SP, (a) as a function of 
oJtt/Ef for Xp = —0.8, (b) as a function of aJn/Ep for Xp = —1.2, and (c) as a function of | Xp | for lOtt = O.lEp. The other 
parameters are the same as in Fig. 3. 



FIG. 5. Phase diagram for X^ = 0.4, Xp = -0.8, Yp = 0.375X^, and = O.AEf- 



FIG. 6. Phase diagram with a cutoff in the particle-particle channel, rj = 0.1. The other parameters are the same as in 
Fig. 5. 
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FIG. 7. Phase diagram with a cutoff in the particle-hole channel, k = 0.1. The other parameters are the same as in Fig. 5. 



FIG. 8. Logarithm of resistance (in arbitrary unit) as a function of temperature (in the unit of Ef) with different cutoffs in 
the particle-particle channel, ?? = 0, 0.02, 0.04, 0.06, 0.08, and 0.1 from the bottom. The parameters are Yi = 0.25, I3 = 0.6, 
and otherwise the same as in Fig. 5. 

FIG. 9. Logarithm of resistance (in arbitrary unit) as a function of temperature (in the unit of Ep) with different cutoffs 
in the particle-hole channel, /t = 0, 0.01, 0.02, 0.03, 0.04, and 0.05 from the top. The parameters are Yi = 0.4, I3 = 0.8, and 
otherwise the same as in Fig. 5. 

FIG. 10. Phase diagram for {Xa^, Xap, Xga, Xsp) = (0.4, -0.4, 0.4, -0.4). 

FIG. 11. Phase diagram for {Xa^, Xap, Xsa, Xbp) = (0.8, -0.8, 0.4, -0.4). 

FIG. 12. Phase diagram for {Xa^, Xap, Xsa, Xsp) = (0.4, -0.4, 0.8, -0.8). 

FIG. 13. Phase diagram for {Xa^, Xap, Xsa, Xsp) = (0.8, -0.8, 0.8, -0.8). 

FIG. 14. Phase diagram for (Xact, Xap, Xsa, Xsp) = (0.5, -0.7, 0.3, -0.9). 

FIG. 15. Logarithm of resistance (in arbitrary unit) as a function of temperature (in the unit of Ef) for different correlation 
strengths. The parameters are {Xau, Xap, Xbu, Xbp) = (0.8, —0.8, 0.4, —0.4), j/i — 0.4, j/3 = 0.3, and w's as before for the 
curve at the bottom. For the other curves, the parameters X's, | y |^'s, and w's are 1.1, 1.2, 1.3, and 1.4 times the above values. 
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